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Type la supernovae (SN la) are generally believed to be the result of the thermonuclear disrup- 
tion of Chandrasekhar-mass carbon-oxygen white dwarfs, mainly because such thermonuclear 
explosions can account for the right amount of 56 Ni, which is needed to explain the light curves 
and the late-time spectra, and the abundances of intermediate-mass nuclei which dominate 
the spectra near maximum light. Because of their enormous brightness and apparent homo- 
geneity SN la have become an important tool to measure cosmological parameters. In this 
article the present understanding of the physics of thermonuclear explosions is reviewed. In 
particular, we focus our attention on subsonic ("deflagration") fronts, i.e. we investigate fronts 
propagating by heat diffusion and convection rather than by compression. Models based upon 
this mode of nuclear burning have been applied very successfully to the SN la problem, and 
are able to reproduce many of their observed features remarkably well. However, the models 
also indicate that SN la may differ considerably from each other, which is of importance if 
they are to be used as standard candles. 



1 Introduction 

Type la supernovae, i.e. stellar explosions which do not have hydrogen in their spectra, but 
intermediate-mass elements, such as silicon, calcium, cobalt, and iron, have recently received 
considerable attention because it appears that they can be used as "standard candles" to mea- 
sure cosmic distances out to billions of light years away from us. Moreover, observations of 
type la supernovae seem to indicate that we are living in a universe that started to accelerate 
its expansion when it was about half its present age. These conclusions rest primarily on phe- 
nomenological models which, however, lack proper theoretical understanding, mainly because 
the explosion process, initiated by thermonuclear fusion of carbon and oxygen into heavier el- 
ements, is difficult to simulate even on supercomputers, for reasons we shall discuss in this 
article. 

The most popular progenitor model for the average type la supernovae is a massive white 
dwarf, consisting of carbon and oxygen, which approaches the Chandrasekhar mass, M^h ~ 



1.39M0, by a yet unknown mechanism, presumably accretion from a companion star, and is 
disrupted by a thermonuclear explosion C3. Arguments in favor of this hypothesis include the 
ability of the models to fit the observed spectra and light curves rather well. 

However, not only is the evolution of massive white dwarfs to explosion very uncertain, 
leaving room for some diversity in the initial conditions (such as the temperature profile at 
ignition), but also the physics of thermonuclear burning in degenerate matter is complex and 
not well understood. The generally accepted scenario is that explosive carbon burning is ignited 
either at the center of the star or off-center in a couple of ignition spots, depending on the details 
of the previous evolution. After ignition, the flame is thought to propagate through the star as 
a subsonic deflagration wave which may or may not change into a detonation at low densities 
(around 10 g/cm 3 ). Numerical models with parameterized velocity of the Imrning front have 
been very successful, the prototype being the W7 model of Nomoto et al. c3. However, these 
models do not solve the problem because attempts to determine the effective flame velocity 
n direct numerical simulations failed and gave velocities far too low for successful explosions 
This has led to some speculations about ways to change the deflagration into a supersonic 
detonation Some of these aspects of supernova models will be discussed in the following 
sections. 
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2 Nuclear burning in degenerate C+O matter 

Owing to the strong temperature dependence of the nuclear reaction rates nuclear burning 
during the explosion is confined to microscopically thin layers that propagate either conductively 
as subsonic deflagrations ("flames") or by shock compression as supersonic detonations at3. 
Both modes are hydro dynamically unstable to spatial perturbations as can be shown by linear 
perturbation analysis. In the nonlinear regime, the burning fronts are either stabilized by forming 
a cellular structure or become fullv turbulent - either way, the total burning rate increases as a 
result of flame surface growth0'@'@. Neither flames nor detonations can be resolved in explosion 
simulations on stellar scales and therefore have to be represented by numerical models. 

When the fuel exceeds a critical temperature T c where burning proceeds nearly instanta- 
neously compared with the fluid motions, a thin reaction zone forms at the interface between 
burned and unburned material. It propagates into the surrounding fuel by one of two mechanisms 
allowed by the Rankine-Hugoniot jump conditions: a deflagration ( "flame" ) or a detonation. 

If the overpressure created by the heat of the burning products is sufficiently high, a hydrody- 
namical shock wave forms that ignites the fuel by compressional heating. Such a self-sustaining 
combustion front that propagates by shock-heating is called a detonation. If, on the other hand, 
the initial overpressure is too weak, the temperature gradient at the fuel-ashes interface steep- 
ens until an equilibrium between heat diffusion and energy generation is reached. The resulting 
combustion front consists of a diffusion zone that heats up the fuel to T c , followed by a thin 
reaction layer where the fuel is consumed and energy is generated. It is called a deflagration 
or simply a flame and moves subsonically with respect to the unburned material 111. Flames, 
unlike detonations, may therefore be strongly affected by turbulent velocity fluctuations of the 
fuel. Only if the unburned material is at rest, a unique laminar flame speed S\ can be found 
which depends on the detailed interaction of burning and diffusion within the flame region El. 
It can be estimated by assuming that in order for burning and diffusion to be in equilibrium, 
the respective time scales, ~ e/w ar id ~ 5 2 / k, where 5 is the flame thickness and k is the 
thermal diffusivity, must be similarll^l: t\ } ~ t&. Defining S\ = 5/t\>, one finds S\ ~ (ra/e) 1 '' 2 , 
where w should be evaluated at T T c (pil). This is only a crude estimate due to the strong 
T-dependence of w. Numerical solutions of the full equations of hydrodynamics including nu- 
clear energy generation and heat diffusion are needed to obtain more accurate values for S\ 
as a function of p and fuel composition. Laminar thermonuclear carbon and oxygen flames at 



high to intermediate densities were investigatednocj, and, using a variety of different techniques 
and nuclear networks, by Timmes & WoosleyEl. For the purpose of SN la explosion modeling, 
one needs to know the laminar flame speed S\ ~ 10 7 . . . 10 4 cm s _1 for p ~ 10 9 . . . 10 7 g cm' 3 , 
the flame thickness 5 = 10~ 4 ... 1 cm (defined here as the width of the thermal pre-heating 
layer ahead of the much thinner reaction front), and the density contrast between burned and 
unburned material p = Ap/p = 0.2... 0.5 (all values quoted here assume a composition of 
Xq = Xo = 0.5 cil) . The thermal expansion parameter p reflects the partial lifting of electron 
degeneracy in thediurning products, and is much lower than the typical value found in chemical, 
ideal gas systems e3. 



3 Hydrodynamic instabilities and turbulence 

The best studied and probably most important hydrodynamical effect for modeling SN la ex- 
plosions is the Rayleigh- Taylor (RT) instability r^a illailg from the buoyancy of hot, burned fluid 
with respect to the dense, unburned material EltE3'Ejl'EJj 13 H 24 l, and after more than five decades 
of experimental aj 
well understoodBt 



mjperical work, the basic phenomenology of nonlinear RT mixing is fairly 
Subject to the RT instability, small surface perturbations grow until 
they form bubbles (or "mushrooms") that begin to float upward while spikes of dense fluid fall 
down. In the nonlinear regime bubbles of various sizes interact and create a foamy RT mix- 
ing layer whose vertical extent h^T grows with time t according to a self-similar growth law, 
/irt = oig{p/2)t 2 , where a is a dimensionless constant (a « 0.05) and g is the background 
gravitational acceleration. 

Secondary instabilities related to the velocity shear along the bubble surfaces^ quickly lead 
to the production of turbulent velocity fluctuations that cascade from the size of the largest 
bubbles ( ~ 10 cm) down to the microscopic Kolmogorov scale, ~ 10~ 4 cm where they are 
dissipated! 24 !! 13 !. Since no computer is capable of resolving this range of scales, one has to resort 
to statistical or scaling approximations of those length scales that are not properly resolved. 
The most prominent scaling relation in turbulence research is Kolmogorov's law for the cascade 
of velocity fluctuations, stating that in the case of isotropy and statistical stationarity, the 
mean velocity v of turbulent eddies with size I scales as v ~ Z 1 / 3 pi). Given the velocity of large 
eddies, e.g. from computer simulations, one can use this relation to extrapolate the eddy velocity 
distribution down to smaller scales under the assumption of isotropic, fully developed turbulence 
cfJ. Knowledge of the eddy velocity as a functio n of length scale is important to classify the 
burning regime of the turbulent combustion front EJJtJ'EJ. The ratio of the laminar flame speed 
and the turbulent velocity on the scale of the flame thickness, K = S\/v{S), plays an important 
role: if A' 3> 1, the laminar flame structure is nearly unaffected by turbulent fluctuations. 
Turbulence does, however, wrinkle and deform the flame on scales I where S\ <C v(l), i.e. above 
the Gibson scale l s defined by Ell S\ = v(l s ). These wrinkles increase the flame surface area 
and therefore the total energy generation rate of the turbulent front 0. In other words, the 
turbulent flame speed, St, defined as the mean overall propagation velocity of the turbulent 
flame front, becomes larger than the laminar speed Si. If the turbulence is sufficiently strong, 
v(L) S> Si, the turbulent flame speed becomes independent of the laminar speed, and therefore 
of the microphvsics of burning and diffusion, and scales only with the velocity of the largest 
turbulent eddyH'B: 



St ~ v(L) . (1) 



Because of the unperturbed laminar flame properties on very small scales, and the wrinkling 
of the flawe on large scales, the burning regime where K ^> 1 is called the corrugated flamelet 
regime BS. 



As the density of the white dwarf material declines and the laminar flamelets become slower 
and thicker it is plausible that at some point turbulence significantly alters the thermal flame 
structure tSBII. This marks the end of the flamelet regime and the beginning of the distributed 
burning, or distributed reaction zone, regimeEl. So far, modeling the distributed burning regime 
in exploding white dwarfs has not been attempted explicitly since neither nuclear burning and 
diffusion nor turbulentmixing can be properly described by simplified prescriptions (see, how- 
ever, Lisewski et al. c3). Phenomenologically, the laminar flame structure is believed to be 
disrupted by turbulence and to form a distribution of reaction zones with various lengths and 
thicknesses. In order to find the critical density for the transition between both regimes, we need 
to formulate a specific criterionJbr flamelet breakdown. A criterion for the transition between 
both regimes is discussed inc3'E3 andEJ: 

Cutoff < £ • (2) 

Inserting the results of0 for S\ and 5 as functions of density, and using a typical turbulence 
velocitv u(10 6 cm) ~ 10 7 cm s _1 , the transition from flamelet to distributed burning can be 
shown c3 to occur at a density of p^is ~ 

10 7 g cm . 



4 Modeling turbulent thermonuclear combustion 

Next we will outline a way by which several of the ideas discussed in the previous sections can 
be incorporated into a numerical scheme to model thermonuclear combustion in SN la. 

Numerical simulations of any kind of turbulent combustion have always been a challenge, 
mainly because of the large range of length scales involved. In type la supernovae, in particular, 
the length scales of relevant physical processes range from 10~ 3 cm for the Kolmogorov-scale 
to several 10 7 cm for typical convective motions. Despite considerable progress in the field of 
modeling turbulent combustion for astrophysical flows (see, e.g., Niemeyer & Hillebrandt@) , the 
correct numerical representation of the thermonuclear deflagration front is still a weakness of 
the simulations. Methods used up to now are based on for the reactive-diffusive flame modelEJ, 
which artificially stretches the burning region over several grid zones to ensure an isotropic flame 
propagation speed. However, the soft transition from fuel to ashes stabilizes the front against 
hydrodynamical instabilities on small length scales, which in turn results in an underestimation 
of the flame surface area and - consequently - of the total energy generation rate. Moreover, 
because nuclear fusion rates depend on temperature nearly exponentially, one cannot use the 
zone-averaged values of the temperature obtained this way to calculate the reaction kinetics. 

Therefore we have decided to use a front tracking method to cure some of these weaknesses. 
The method is based on the so-called level set technique which was originally introduced by 
Osher and Sethianca They used the zero level set of a n-dimensional scalar function to represent 
(n — l)-dimensional front geometries. Equations for the time evolution of such a level set which 
is passively advected by a flow field are given in Sussman et al. e3. The method has been 



extended to a 
detonations™' 



fhe tracking of fronts propagating normal to themselves, e.g. deflagrations and 
In contrast to the artificial broadening of the flame in the reaction-diffusion- 
approach, this algorithm is able to treat the front as an exact hydrodynamical discontinuity. 
We will demonstrate that such a method can be applied to the supernova problem and, in 
addition, we will show that even if one attempts to model the physics of thermonuclear burning 
on unresolved length scales well by physically motivated "Large Eddy Simulations" (LES), one 
still has to perform calculations with very high spatial resolution in order to verify the correctness 
of the employed model for the sub-grid scales. 
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Figure 1: Snapshots of the temperature and the front geometry at 0.4s for model B5 of Reinecke et al. (left figure) 
and turbulent velocity fluctuations on the grid scale (right panel). The position of of the front is indicated by the 

dotted curve. 
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Figure 2: Snapshots of the temperature and the front geometry at 1.05s, taken for the low-resolution model (left 

figure) and high-resolution run, respectively. 

5 Application to the supernova problem 

We have carried out two-dimensional numerical simulations in cylindrical rather than in spherical 
coordinates, mainly because it is much simpler to implement the level set on a Cartesian (r,z) 
grid. Moreover, the CFL condition is somewhat relaxed in comparison to spherical coordinates. 
The grid we used in most of our simulations maps the white dwarf onto 256x256 mesh points, 
equally spaced for the innermost 226x226 zones by A=1.5-10 6 cm, but increasing by 10% from 
zone to zone in the outer parts. The white dwarf, constructed in hydrostatic equilibrium for a 
realistic equation of state, has a central density of 2.9-10 9 g/cm 3 , a radius of 1.5-10 8 cm, and a 
mass of 2.8-10 33 g, identical to the one used by Niemeyer and Hillebrandt c3. The initial mass 
fractions of C and O are chosen to be equal, and the total binding energy turns out to be 
5.4-10 50 erg. At low densities (p < 10 7 g/cm 3 ), the burning velocity of the front is set equal 
to zero because the flame enters the distributed regime and our physical model is no longer 
valid. However, since in reality some matter may still burn the energy release obtained in the 
simulations is probably somewhat too low. 

Here we present only the results for two models, one in which nuclear burning was ignited 
off-center in a blob and a second onein which initially five blobs were burned as an initial 
condition, and refer to Reinecke et al. ^3 for simulations with other initial conditions. 



First, in Fig.p] a snapshot of the "five-blob" model B5 is shown at t = 0.4s. The left panel 
gives the position of the burning front, the temperatures (in gray-shading), and the expansion 
velocities. The right panel shows the distribution of turbulent velocity fluctuations. We find 
that, in accord with one's intuition, most of the turbulence is generated in a very thin layer near 
the front. Since in the limit of high turbulence intensity the nuclear flames propagate with the 
turbulent velocity it is obvious that this propagation velocity exceeds the laminar flame speed. 
However, for most of the initial conditions we have investigated this increase was not sufficient 
to unbind the star. In fact, model B5 of Reinecke et al. H was the only one of the set that did 
explode. 

Moreover, we show the results obtained with 3 times higher resolution for comparison. Figj2| 
gives a snapshot taken at 1.05s after ignition for the low-resolution run (left figure) and the 
high-resolution run, respectively. Although the increase in spatial resolution is only a factor 
of 3, the right panel shows clearly more structure. This is an important effect because in the 
flamelet regime the rate of fuel consumption, in first order, increases proportional to surface 
area of the burning front. The net effect is that the low-resolution model stays bound at the 
end of the computations, whereas the better resolved model explodes with an explosion energy 
of about 2xl0 50 erg. Fig.^ also demonstrates that the level-set prescription allows to resolve the 
structure of the burning front down almost to the grid scale, thus avoiding artificial smearing of 
the front which is an inherent problem of front-capturing schemes. We want to stress that this 
gain of accuracy is not obtained at the expense of smaller CFL time steps because in our hybrid 
scheme the hydrodynamics is still done with cell-averaged quantities. 

Unfortunately, the noticeable increase in the total energy release for the better resolved 
simulation is a strong indication for problems with our sub-grid scale model. Ideally the increased 
flame surface should have been exactly balanced by a smaller St, and the total amount of 
released energy should not have changed much. This aspect has to be investigated further 
before convincing predictions about the explosion energetics can be made. Nevertheless the 
results of the better resolved model are expected to be more accurate since a larger scale range 
is covered by direct simulation and the sub-grid model should be less important. 

Recently wehaveper formed the first three-dimensional supernova simulation based on the 
presented numerical models. The calculation was carried out on a cartesian grid of 256 3 cells 
with the same spacing as described for the 2D case. To save computation time, only an octant of 
the star was simulated and reflecting boundaries were imposed on the coordinate planes, which 
can be interpreted as an eightfold symmetry. This time the star was ignited in the center, with 
a perturbation of the radius of the ignited region depending on This situation is the three- 
dimensional equivalent to the model C3 presented inE3. Figure || shows a snapshot of the front 
geometry after 0.5s. It is evident that the initial axial symmetry has disappeared completely, 
resulting in a larger flame surface than in the 2D case where the symmetry was artificially 
enforced. As a consequence, the energy release is considerably higher in three dimensions and is 
more than sufficient to unbind the star, whereas the corresponding 2D model was only marginally 
unbound. 

We want to stress again, however, that the numerical simulations described here, although 
leaving behind unbound white dwarfs in several cases, would not fit the observed spectra of SN 
la, mainly because the expansion velocity of the partially burnt gas is too low. 

6 Summary and conclusions 

In this article we have outlined our present understanding of the explosin mechanism of Type 
la supernovae. From the tremendous amount of work carried out over the last couple of years 
it has become obvious that the physics of SNe la is very complex, ranging from the possibility 
of very different progenitors to the complexity of the physics leading to the explosion and the 



Figure 3: Flame surface in the three-dimensional C3 model at t=0.5s. The integral scale of the flame is about 

1.5 ■ 10 8 cm. 

complicated processes which couple the interior physics to observable quantities. None of these 
problems is fully understood yet, but what one is tempted to state is that, from a theorist's 
point of view, it appears to be a miracle that all the complexity seems to average out in a 
mysterious way to make the class so homogeneous. In contrast, as it stands, a safe prediction 
from theory seems to be that SNe la should get more divers with increasing observed sample 
sizes. If, however, homogeneity would continue to hold this would certainly add support to the 
Chandrasekhar-mass scenario discussed in this article. On the other hand, even an increasing 
diversity would not rule out Chandrasekhar-mass progenitors for most of them. In contrast, there 
are ways to explain how the diversity is absorbed in a one parameter family of transformations, 
such as the Phillips relation H'i or modifications of itSS. 

As far as the explosion/combustion physics and the numerical simulations are concerned 
significant recent progress has made the models more realistic (and reliable). Thanks to ever 
increasing computer resources 3-dimensional simulations have become feasible which treat the 
full star with good spatial resolution and realistic input physics. Already the results of 2- 
dimensional simulations indicate that pure deflagrations waves in Chandrasekhar-mass C+O 
white dwarfs can lead to explosions, and going to three dimensions, because of the increasing 
surface area of the nuclear flames, further adds to the explosion energy. On the side of the 
combustion physics, the burning in the distributed regime at low densities needs to be explored 
further, but it is not clear anymore whether a transition from a deflagration to a detonation in 
that regime is needed for successful models. In fact, according to recent studies such a transition 
appears to be rather unlikely. 
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